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Abstract 

The correction to the wave function of the ground state in a hydrogen-hke atom due to an 
external homogenous magnetic field is found exactly in the parameter Za. The j = 1/2 projection 
of the correction to the wave function of the nsi/2 state due to the external homogeneous magnetic 
field is found for arbitrary n. The j = 3/2 projection of the correction to the wave function of the 
ns 1/2 state due to the nuclear magnetic moment is also found for arbitrary n. Using these results, 
we have calculated the shielding of the nuclear magnetic moment by the nsi/2 electron. 
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I. INTRODUCTION AND CONCLUSIONS 



Nowadays, the measurements of magnetic properties of a bound electron or nucleus in sim- 
ple atoms have reached an impressive level of precision. This progress stimulates theoretical 
investigations of subtle effects which contribute to the magnetic moment of a hydrogen-like 
atom. To some extent, a nucleus can be considered as a source of the Coulomb field dis- 
torted at small distances due to the nuclear size effect. However, the present experimental 
accuracy leads to the necessity to go beyond this approximation. For a spinless nucleus one 
has to take into account the recoil effects, which lead to the corrections in mg/mTv (w^e is 
the mass of electron and m^r is the nuclear mass). For the nucleus with a nonzero spin and 
magnetic moment, one also has to take into account the interaction of the nuclear magnetic 
moment with the electron and external magnetic field B. E.g., for a spin-1/2 nucleus, these 
effects can be described by the Hamiltonian [l| 

6 6 

T^nij = 7^ — g'e 3e - B - - — Sat ■ S + A ■ Sat + higher order terms , (1) 
zme Irrtp 

which is defined in the space of spin and angular moment variables for the hyperfine compo- 
nents of the nlj state in the hydrogen-like atom. Here, e is the elementary charge (positive), 
rrte/p is the electron/proton mass, is the electron total angular momentum, is the 
nuclear spin and A is the hyperfine-interaction parameter. The relativistic units in which 
h = c = 1 are applied through the paper and = a ~ 1/137 is the fine structure constant. 
The bound g factors of the electron and nucleus (g'^^/j^) are somewhat different from their 
free values gf^j^ due to bound effects which increase with Z. 

Study of such a bound correction for the nuclear g factor (i.e. the screening or shielding 
effect) for certain states was done some time ago in . This result was checked for the 
Is state in Ml. An important detail of the latter calculation is the related correction to 
the wave function which was partly obtained for the Is state in 3, M for other purposes. 

I nri 

Actually, one of the purposes of was to verify the previously obtained [5|, la] results for 
the wave function. 

In this paper, we consider the effect of shielding of the nuclear magnetic moment by the 
electron in nsi/2 state, which results in the difference between g'^ and g^^\ For this purpose, 
we calculate the corrections to the electron wave function caused by both an external mag- 
netic field and the nuclear magnetic moment. The operators corresponding to these effects, 
have vector character and the correction to the wave function of an s state is a mixture of 
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j = 1/2 and j = 3/2 parts. The form of these parts is quite different. In particular, the 
j = 1/2 projection of the correction due to the external magnetic field for the ns state has 
the form of some polynomial operator acting on the unperturbed wave function, Eq. (1151) . 
The j = 3/2 projection for this correction can not be represented in the similar form. For 
the specific case of Is state, this correction is calculated below. For the correction due to 
the nuclear magnetic moment the situation is the opposite: the j = 3/2 projection has the 
form of some polynomial operator acting on the unperturbed wave function (see Eq. f|T6|) , 
while the j = 1/2 projection can not be expressed in such a simple form (see Appendix). 

Using the derived corrections, we have calculated the shielding of the nuclear magnetic 
moment in nsi/2 state. The expansion of our result (see Eq. f[T7|) for a result exact in Za) is 

a(Za) 



q'n = 9n^ X 



3^2 



132n-35,^ ,2 370^3 + 342^2 - 492n + 69 , ^ 




1320872^ + 15048n^ - 8552^3 - 25320^2 + 17064ra - 1641 , ^ 

H 7. (Za) + . . . 

1728^6 ^ ' 

For the specific case of the ground state, this expression agrees with those in [2l, ISj. Our 

method of calculations is, inpart, similar to theirs and the related results on the wave 

function in our paper and in are also consistent. 

The paper is organized as follows. In Section [TTl we calculate the correction to the wave 

function of the ground state in a hydrogen-like atom due to an external constant magnetic 

field. We check this correction by calculating the magnetic polarizability and the shielding 

of the nuclear magnetic moment in the ground state, which are known. In Section IIIII we 

calculate the shielding of the nuclear magnetic moment by the nsi/2 electron. 

II. CORRECTIONS TO THE WAVE FUNCTION OF THE GROUND STATE IN 
THE EXTERNAL MAGNETIC FIELD 

In the present Section we calculate the correction to the wave function of the ground 
state due to an external constant magnetic field B. 

The unperturbed wave functions of the ground state with = ±1/2 have the form 



il>± = C [mry ^ exp [—Zamr) 



Vl + 7 V^i 



VI - 7 i{cr ■n) Lp± 
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where 



n — r/r , C 



V87rr(27 + 1) ' 











: 

















Let the external magnetic field be directed along ;2-axis. The Hamiltonian has the form 

if = a • (p - eA) - Za/r + 7°m ^ Hq- e{ct ■ A) , 

A^- B xr. 
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The states ip^ form a diagonal basis for the perturbation. The correction to the wave 
function -0+ is expressed as 

eB 

Sijj — —eGr (a • A) ■0+ = —Gr (cc x r)^ ■0+ , 
where Gr is the reduced Green function determined as 

~l'fa)<<*±l| (2) 



and £ = m7 is the energy of the ground state. In order to determine b'll) explicitly, we first 
separate parts with j = 1/2 and j = 3/2. They can be expressed in terms of the four 
functions gi, hi, g2, /i2 yet to be determined: 

Sips'^ = C - — {mry^^ exp {—Zamr) [ ^ ^ 9i{ )v+ i 



2m 



VI - 7 ^1 (^) i{cr-n) 



^'4'^^ — C -—- {mry ^exp(— Zamr) [ 7 5'2( ) [ ( ) ^ I (4) 

2 \^ Vr^/i2(r)z[^-(o--r)/3](^+ 

Using the fact that Hq commutes with j and that the projection operators on the sub- 
spaces with j = 1/2 and j = 3/2 in our case can be represented as 

Pv.^j-f ^3/. = £-i, (5) 
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we have 



Sipg"^ = —C - — Gr {mry exp {— Z amr) i ^"^^ | , (6) 

"^"^ y VI +7^ (cr ■ n) V9+ 



^i^B — ^ — {mry ^ exp {—Zamr) ( 1 . (7) 

2 \^ Vr+7^ [2: - (o- ■ r) /3] J 

Acting by the operator e — Hq on both sides of Eqs. (jHl ED and using Eqs. ([31 S]), we 
obtain two systems of equations corresponding to j = 1/2 and j = 3/2. 



A. The j = 1/2 projection of the correction to the wave function of the ground 
state 

Taking into account that 

(£-Ho)G, = l-|^+)(^+|-|^-)(^-| , 
we obtain the following system of equations for the contribution of j = 1/2: 

V ^ J 3(1+7) 3 ' 

w A + 7 ,^ Za(27 + 1) 2mr 

The boundary conditions can be easily determined from the representation of the reduced 
Green function in terms of the eigenfunctions of the Hamiltonian. Since the wave functions 
of the nsi/2 states behave as r'^^^ at the origin, we conclude that the functions gi and hi 
should be regular at this point. 

Let us pass to the functions gi± = gi± hi. For gi_ we obtain a closed equation: 

2 27 + 1 Amr 



g'i_ + 2 (^^ - mZa^ gi. 



3 Za 



The general solution of this equation is 



2r 

gi- = \- Ci{mr) exp [2 Zamr) . 

3Za 

Due to the boundary conditions, Ci = 0. Using Eq. ([H]), we obtain 

Za (27 + 3) 
3(1+7) ^ + 

hi^-^^^P^r + C.. 
3(1-7) 
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The constant C2 is completely determined by the condition of the orthogonality of Sip'^/'^ 
and Finally, we obtain 

Za (27 + 3) 27 + 1 



9i 



3(1 + 7) 3m 



Za (27 - 1) 27+1 

hi = — r , 

3(1 — 7) 3m 

which agrees with 

B. The j = 3/2 projection of the correction to the wave function of the ground 
state 

We obtain the following system of equations for the contribution of j = 3/2: 
92 = — \ mZa j 92 — \ h rnZa j 112 + mr , 



h'2 = ^— — mZa^ g2 + ^ + mZa^ h2 — 



2 — mr . 



Again, moving to g2± = (72 ± ^2, we have 

92- = -rg2+ , (9) 
3 

r (^2+ + (27 + 1 — 2Zamr) (^2+ 92+ = — 2mr . 

The first equation determines the function (72- via g2+. The general solution of the second 
equation is 

c,2+ = Ciimr)-"'^ iFi (-z/ - 7; 1 - 2z/; 2Zar) + C2 {mrY'^' iFi (z/ - 7; 1 + 2u; 2Zamr) 

\2 



2 ( TTIT 1 

^ 2i^2 (2, 1; 3 + 7 + i^, 3 + 7 - u; 2Zamr) 

47 + 1 



where i/ = y4 — (Za)'^. 

Again, Ci = due to the boundary conditions at the origin. The constant C2 is fixed 
to suppress the exponential growth of (72+ as r — *• cxo. In order to determine it, we use the 
following integral representation of 2F2 (see, e.g. [Sj): 

2F2(2,l;a,6;x) 



(a-l)(a-2) (6-1) 



[ dzz{l-zy-^ [ dtt^-^exp[z{l-t)x] . 
Jo Jo 
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We rewrite the integral over t as dt = dt — dt. The first integral results in 
the F-function, and we obtain the integral representation for the iFi function. Making the 
substitution t — >• 1/t in the second term, we arrive at 

(2, 1; a, b; x) = ^i-^ r (a) F W T (3 - fe) (3 - 6; 1 + a - 6; x) 

T [1 + a — b) 



(a-l)(a-2)(6-l) / dzz{l-z 



\a—3 







[ dtt'^exp [z{l - l/t)x] 
Jo 



In our case, a = 3 + 7 + z/, 6 = 3 + 7 — z/, x = '2Zamr and the first term contains 
i-Fi (z/ — 7; 1 + 2z/; 2Zamr). Note that the second term in the right-hand side does not grow 
exponentially as r ^ 00. Thus, we have 

^ _ 2r(3 + 7 + ^)r(3 + 7-^)r(z/-7) 

m(l + 47)r(l + 2z.) ^^^"^ 

and the following integral representation for g2+ is valid: 

g2+ = 2mr^ (1 + 7 + i/) /" dz z {1 - zf^" f dt t-^^"-^ exp (1 - 1 /t) 2Zamr] . 

Jo Jo 

The function g2- is determined by Eq. iQ. Note that, in contrast to the j = 1/2 
projection, this correction cannot be expressed via elementary functions. This circumstance 
complicates the derivation of the j = 3/2 projection of the correction due to the external 
magnetic field in ns states with arbitrary n. 



C. Magnetic polarizability and nuclear magnetic moment shielding in the ground 
state of the hydrogen-like atom 

As an application of the obtained correction to the wave function, let us first calculate 
the magnetic polarizability of the hydrogen-like atom, which is related to a correction to the 
energy level quadratic in a magnetic field. 

Though the Hamiltonian ([1]) is linear in the magnetic field, its eigenvalues are not. One 
can consider them at the low-field limit (i.e. the field weak in comparison to the hyperfine 
splitting), which contains terms non-linear in B (see, e.g. 

Here we derive the polarizability for another case, when the hyperfine splitting can be 
neglected (that is a kind of 'medium' field, which is weak comparing to the gross structure 
splittings, but strong in respect to the nuclear hyperfine field). Once we neglect the hyperfine 
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splitting, the magnetic polarizability (3 can be presented in the following form: 
8E=^= ' j d'r^i (a X r), Gr (a x r), ^+ 

Using Eqs. ©-(I?]), we obtain 

(3 = f {mrf^-^ exp{-2Zamr) {g2+ - 3gi+) . 

9m^ J 

Taking the integrals, we have 



^_ «(27 + l 



(z/ + 7) (27 + 3) (z/ + 7 + 1) , 

I ^ /M 1^ M ^I^ ^3ir2(2,2,z/-7;z/ + 7 + 4,z/-7 + l;l) + 47-2 



36m3 (1-7) [ 3 (i/ + 7 + 2) (z/ + 7 + 3) 

(10) 

This result agrees with Qj. 
In the limit Za <^ 1 we have 

(3 " ^ ^ 



m3 \2{Zay 3 ' 

Let us now calculate the shielding of the nuclear magnetic moment in the ground state. 
The correction to the energy of the ground state in the external constant magnetic field 
reads 

d^r tpl — ^Gr {a x r)^ ip+ , 

where we introduce the shielding constant k as 

g' = {l + k)g^'^ . 
Again, taking the integral, we obtain 



in agreement with 



AaZa (I 1 2 3 ■ 

9 V3 6(1 + 7)^7 27-1,'' ^ ' 
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III. SHIELDING OF THE NUCLEAR MAGNETIC MOMENT IN nsi/2 STATES 



Let us now calculate the corrections to the wave functions of the nsi/2 state for an 
arbitrary n. Using the projection operators ([S]), one can represent the correction to the 
energy in the external constant magnetic field as 

= /inuci5 j dh ^^i^^^Gr (P3/2 + P1/2) (« X r)^ 

= 2/i„,cie j d\ijl^^^fi^6^]l' + Be j dh {5ipy^)\cx x r) _^,p+ . (12) 

Here 5ip^ and 5ipB denote the first-order corrections to the wave function due to the 
magnetic moment of the nucleus and due to the external field, respectively. As we shall 

3/2 1/2 

see, this representation is convenient for the calculations, since both 5ip^ and Sil)^ are 
expressed via the action of some operators, polynomial in r, on the wave function of 
the nsi/2 state. This wave function has the form 



, a\/m + £09+ 

V'n.,/, = C {mrV exp (-Ar) | | , (13) 

ih\/m — e {cr ■ n) if^ 

where C is a normalization constant and the polynomials a and b satisfy the following system 
of equations 

a' = (-A - ^ (m - £)) 6 + (a + a , (14) 

6'=(-A + ^(m + .))a+(A-l±l)6, 

a' =h-^)a,-(l + ^-2X\a. 
\r Xr J \r Ar 



7 Zae^ _|_ ^''^ _|_ ^^^^ 



0+ = I \- — ; — I 0+ + I - H ; — I o- 

r Ar 

Here e is the energy of the state, A = \Jm? — e^, a± = a ± 6. For lsi/2 state, a = 6 = 1. 
The idea of our derivation of the corrections to the wave functions is to use these differential 
equations rather than the explicit form of the functions a and h. 
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1/2 

A. Calculation of 

1 /2 

We start from the derivation of the correction to the wave function 5ip^ due to the 
presence of the magnetic field. Again, using the projectors ([5]), we obtain 



1/2 eB I h\Jra — e </3+ \ eB 

^i^B = -C iz—Gr {mry exp {-Zamr) = —Gr {-ij ■ r) i^ns^^ ■ 

\ a\/m + ei {(T ■ n) Lp+ • 

Using the following identities 

3 i 
-im'j ■ r = - + ■ L + -[oL ■ r, Hq] 



/2 



= (1-7°) ^ns, 

and the fact that Gripns^/2 = 0, we arrive at 



The first term in the braces is conveniently represented as 
^ , ^ , 1/^3 d 

Using Eqs. (fT3l) and ([HD, we have 



c- ,1/2 eB f 1 7° / £ 1 Za\ . ■ \ , /-,^\ 

3m V 2m m Vm 2 mr ' ' ' 



B. Calculation of 6ip'^^'^ 

Let us find the correction to the wave function due to the magnetic field induced by the 
nuclear magnetic moment. 

by^m — e [z (cr ■ n) — r/3] 



6ijl/' = C e/i„,ei Gr {mry-^ exp (-Ar) 



a-y/m + e i[z — {cr ■ r) /^] tp^ 



C e/inuci {mry ^ exp (-Ar) 



g\/m + e [2; (cr ■ n) — r/3] 



/ii/m — e i[z — {cr ■ r) /?)] 
Acting by the operator e — Hq, we obtain the following system 



7 + 2^ Zam 

Ar ^ 'J ' r^ 
Zam \ /7 - 2 \ b 
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Adding and subtracting these equations, we have ior g± — g ± h 

1 / Zam 
r V A 



9'^ 



2 



1 / Za \ a_ 



1 / Zam 



X 



+ 2 U+ 



A 



0\ 7 ^« \ «H 



and, after ehminating of g_, 



-rg'l - (1 + 27 - 2Ar) - (^2X{n -l)-^^g^^ 



a_ + a_| 



Keeping in mind that a± are some polynomials, we may search for the solution of this 
system in the form of power series. It is obvious, that the number of terms in these 
power series, if finite, grows with n. Instead, we search for the solution in the form 
g+ — a+ J2k ^kf^ + ^- X^fe ^kf^ compare the coefficients at r^a±. 

We arrive at the following system of recurrence relations 

-{k^ + 2k (n-l + 7)-3)Cfc-2A; {l - ^ Sk,-i - 2 {k - 1) XC^-i , 

(Z \ 
1 + j Ck = Sk,-i + 2kXDk-i . 

For A; = — 1 we have a system 

2(n + T)C_, + 2('l-^)D_, = l 



X ^ / Zam\ _ 

2 (2 - n - 7) L>_i + 2 M + j C_i = 1 



with the solution 



For A; = 0, we have 



m 



2XZa 



C - ^^C 



and, obviously, the coefficients for higher k are equal to zero. 
Thus, we obtain 

m — e 



9+ 

9- 



2Za 
m — e 
2Za 



1 2\ 1 



A7"+-li-3'"- 
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and 



9 
h 



m — e 

Za 
m — e . 
3Za 



1 . 1 ' 



■b. 



In matrix notations the result can be presented as 



^/^nucl 



m 
3Za 



r 3 



1 - 



m 2mr 



;i + 7o) h ■ n 



(16) 



For the particular case n = 2 the expression agrees with Eqs. IA8I and IA9I derived indepen- 
dently. 



Total result 



Using Eqs. (|T2l) . ( ITSj) and ( |T6l) . we obtain the following result for the shielding coefficient: 



2a , I 1 
k = — — {nsi/2 \ — 7 
9 mr 



2i 



7 • r 



m 



e 1 
-?7 ■ r 



7 



3Za; 2Za 

The first two terms correspond to the contribution of j = 1/2 while the last two terms to 
that of j = 3/2. This matrix element can be evaluated using the recurrence relations given 
in [31. The final result reads 



k = - 



2a 
2a 



A 

Za V m 



A 



+ 



4 (2e + m) /A 
rwy (47^ — 1) \m 



(m — e) {m + 2e) ^ m — e 



QZam? 



2Zam 



m J m7 (47^ — 1) 
and its expansion for small Za is 



4 {2e + m) ^{m- e) (5m + As) 



SZarv? 



(17) 



a{Za) 
3n'^ 



1 + 



132n - 35 



36n^ 



(Zay + 



370^3 + 342^2 - 492n + 69 
72^1^ 



(ZaY 



13208^5 + 15048^^-8552^3- 25320722 + 17064n- 1641 ^ ,0, 
+ . (^a)^ + O ((^a)^) 



(18) 



1728n6 

Only the first term of this expansion has been known for arbitrary n, while the other terms 
were known only for the Is state. 

Note that the found correction corresponds to the following term in the effective Hamil- 
tonian 

6H = -k (-B^„^^i) . 



12 



Another structure, linear in the nuclear magnetic moment and magnetic field 

(^Kucl X -^c]) 

is T-odd and does not appear, which is consistent with the results |6|] (cf. also 

0,3). 
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APPENDICES 

APPENDIX A: DIRECT CALCULATION OF 8^)]!'^ FOR 2s STATE 



The wave function correction due to the nuclear magnetic moment can be calculated 
using the method developed in [7]. For the j = 1/2 part of the correction to the ground 
state it was obtained in [5] confirming the result of j^. For the 2s state one can obtain^ 

-i(cr • n) Y^^r) f^i/2,o,»n(n) 



1/2 



where 



2s 



-i{(T ■ n) g2s{r) ^]l/2,o,i,(r^) 
is the unperturbed state, 17^ ; ^(n) are the spherical spinors. 



X 



is) 



1 r/2A3 

1 — 47^ \m'^'j 
2Zam{2e/m + 1) 



is) 



2s 



3 
r 

F 

3 



2s 



-3m{e/m + 1] 



2Za\ 

92s 

r ) 



2A3 

m"^7 r 



1-472 
2Zam(2e/m + 1 



92s + ?>m{e/m - 1) 



r J 



f'. 



2s 



-G 



(Al) 



(A2) 



(A3) 



(A4) 



^ Note a different representation of the wave function correction as compared with the above text. Here we 
follow [a, mi, where the related results were obtained for the Is state in a similar manner. 
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{Zamy - 2{Zamf\ + ZamX^ - 2A^ 2Zam ^{l + 27) 
2(Z«m- A)A3 ^ A 

2(Zam)^ + A{Zamf\ - b{Zamf\^ + SZamA^ - 2A^ 4A ^(1 + 27) 



F = A^(2Ar)^e~^V2A/m + Za 
+ 
4 



ZamA(Za;m — A)^ 



Zam — A 



Ar 



(Ar)2 ^ — log(2Ar) + Arlog(2Ar) 



Zam — A A 
G = -N{2Xrye-^''^2X/m - Za 



Zam — A 

{Zamf + 5{ZamYX^ - SZamX^ - AX^ 
2Zam{Zam — X)X^ 



(A5) 



+ 



2(Zam + 2A) V(l + 27) 
A 



^ -2{Zamy + 2{ZamfX- {ZamfX"^ -SZamX^ + 2X^ 4AV(l + 27)' ^ 



+ 



2A 



Zam — X 



(XrY - 



N 



ZamX{Zam — A)^ 
2 2(Zam + 2A) 



Zam — A 



A 



log(2Ar) + 



4A 



Zam — X 



Ar log(2Ar) 



(A6) 



A=^ 



1 + 27 



2((Zam)2 - 2A2) V Zam{Zam + A)r(l + 27) ' 
^{z) = T'{z)/T{z) . 

The result for the j = 1/2 part of the 2s shielding obtained using this expression agrees 
with Eq. (fTTl). 

The same method can be used for the calculation of the j = 3/2 part of the correction 
for any ns state and it yields 



-i{a- ■ n) Fi?(r) ^3/2,2,™!^) 



(A7) 



where 



Za[ 3 ^"^ + r 



1 e — m 
Za 3 



(A8) 
(A9) 



which agrees with Eq. ( |T6|) . 
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